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We analyse the classical model of competition between three species studied by 
May and Leonard (SIAM J Appl Math 29 (1975) 243-256) with the approaches of sin- 
gularity analysis and symmetry analysis to identify values of the parameters for which 
the system is integrable. We observe some striking relations between critical values 
arising from the approach of dynamical systems and the singularity and symmetry 
analyses. 

1 Introduction 

In a classic study of a model of competition among three species May and Leonard ^5] 
demonstrated the dramatic change in the qualitative behaviour of the model in simply 
going from two species to three species. The Gausse-Lotka-Volterra [7j HI] model for 
competition among n species is 



where Ni (t) is the size of population i at time t, r\ is its intrinsic growth rate, a y - the 
coefficient representing the effect on its growth rate due to species j and overdot denotes 
differentiation with respect to time. May and Leonard restrict the number of competing 
species to three and make some assumptions about the parameters in the system to reduce 
the system to one which is susceptible to analytic treatment in the main. The critical point 
in choosing n = 3 is not the smallness of the number, but the potential for dramatic change 
in the behaviour of the system in going from n = 2 to n = 3. When n = 2, the autonomous 
system, with n = 2, can be reduced to a single first-order equation and is integrable. 
This is not automatically the case for n = 3. Indeed the potential, if not its realisation, for 
chaos exists. The parameters, rj, i = l,n, are taken to be equal and then set at unity by 
a rescaling of time. Rescaling of the independent variables enables the diagonal elements 
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of the quadratic terms to be set at unity. Finally the interaction coefficients are limited 
to just two by the assumptions that species i + l(mod3) affects species i, i = 1,2,3, and 
species i + 2(mod3) affects species i, i = 1,2,3, in the same way. The model system is 
then 

x = x (1 — x — ay — (3z) , 



z = z (1 — ax — (3y — z) . 

We remark that the simplifications made to the values of the parameters are not as restric- 
tive as one would imagine. For example in grasslands the reproductive rates of different 
species of ungulants of similar size are expected to be similar and the coefficients of com- 
petition likewise. Indeed under good grazing conditions the ciy would be anticipated to 
be low and of comparable magnitude. The mathematical attraction of this model is that 
the community matrix, videlicet 



is a circulant matrix for which an explicit formula for the eigenvalues exists With the 
entries of A as indicated its eigenvalues are 



May and Leonard give the equilibrium points of system 1)1. 2|) as (0, 0, 0); (1, 0, 0), (0, 1, 0) 
and (0, 0, 1); (1 - a, 1 - 0, 0) fa, (1-0,0,1- a) fa and (0, 1 - a, 1 - (3) fa, where 7 = 
1 — a/3; (1, 1, 1) / (1 + a + 0) for zero-, one-, two- and three-population equilibria. In this 
paper we investigate the properties of system ((1.2(1 from the approaches of singularity 
analysis and symmetry analysis. We emphasise that the thrust of our investigations is 
the integrability of system p. 2(1 and not its qualitative behaviour for which the methods 
of dynamical systems are well-suited. The singularity analysis is directed towards the 
determination of the existence of solutions which are analytic. Symmetry analysis leads 
towards invariance of the system under infinitesimal transformation so that in the presence 
of a suitable number of symmetries the solution of the system may be reduced to a sequence 
of quadratures or the existence of three functionally independent invariants from which 
the solution follows by a process of elimination of variables. In the case of the latter the 
elimination may be only local through the use of the Implicit Function Theorem. Equally 
the performance of the quadratures may not be possible in closed form or lead to analytic 



Before we begin any analysis we observe that under a constraint upon the parameters 
a and system 1(1.2(1 is an example of a decomposed system since, if we add the three 
equations, we have 




(1.2) 



A 



1 a 
1 a 
a (3 1 



(1.3) 



Aj = l + a + 0, 
A 2 ± = - 2-a- 0±iV3(a- 0) 



solutions. 



(x + y + z)- = (x + y + z) - {x 2 + y 2 + z 2 + (a + (3){xy + yz + zx)} 



(1.4) 
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Clearly the constraint a + (3 = 2 enables us to write (|1.4|) as the composed system 



u — u 



where x + y + z = u, which is readily integrated to give the invariant 
1 



h 



1 



x + y + z 

and this can be rearranged as the analytic solution 

x + y + z 



(1.5) 



(1.6) 



(1.7) 



for the total population. (Recall that time was rescaled; this explains the simplicity of 
the time dependence in (|1.7|1 ). One of the attractive features of decomposible systems is 
that the composed equation, particularly in the case of systems of first-order differential 
equation of the type usually encountered in modelling, is usually integrable so that an 
invariant exists and the dimension of the system is effectively reduced by one |H1 E3 • 



2 Singularity analysis of system ( 11.21 ) 

We follow the standard method of singularity analysis 2 and determine the leading-order 
behaviour by setting x = At p , y = Br q and z = Ct t , where r = t — in and to is the 
location of the putative movable pole, in system (|l,2j) to obtain 



pArP- 1 



rCr 



r-1 



At p (1 - At p - aBr q - fiCr r 
Br q (l- (3At p - Br q - aCr r 
Ct t (1 - aAr p - (3Br q - Cr r 



from which it is evident that the linear terms of the right hand side are not to be con- 
sidered for the determination of the leading-order behaviour or of the resonances. On the 
assumption that the leading-order behaviour assumed does in fact represent polelike be- 
haviour in the three dependent variables the requirement of balance of the terms reduces 
to just —l,p,q,r from which it is evident that p = q = r = —1. The possibility that 
one or other of the exponents differs from —1 cannot be entertained since the exponent 
would then be nonnegative and the singularity analysis does not admit such a possibility 
for integral leading-order behaviour. One could imagine the introduction of branch point 
singularities with fractional exponents, but this leads us away from the standard analysis. 

With the common exponent of the leading-order behaviour being —1 the coefficients of 
the leading-order terms are found from the solution of the system 



1 

(3 
a 



a 
1 



(3 
a 
1 





A 




1 




A 




B 




1 




B 




C 




1 




C 



1 



1 + a + P 



2 The reader is referred to Ramani et al [22\ and Tabor |2.'j| for an account of the details of the application 
of the Painleve Test and implementation of the ARS algorithm. 
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and we note that the solution is the same as the location of the interior equilibrium point 
of the system (|l,'2j) . 

The resonances are determined from system (jl.2|) (with the linear term omitted) by 
the substitutions 



At' 1 + MT r ~\ 



y 



Bt^+Nt 1 - 1 . 



Ct- 1 + St 



r-1 



where r denotes, as usual, the resonance and should not been confused with the usage 
above as one of the exponents of the leading-order behaviour. The terms linear in M, N 
and S give the eigenvalue problem 



(2.1) 



r + A 


aA 


pA 




' M ' 




' " 


(3B 


r + B 


aB 




N 







aC 


pc 


r + C 




S 








from which it follows that 



-1 and r 2 ± = - 



a + P -2±iV3(a- f3) 
1 + a + P 



(2.2) 



If we look at the eigenvalues of the Jacobian matrix of system ()1.2|) at the interior 
equilibrium, we find that 



A, 



(1 + a + P) n, 



1,2±. 



May and Leonard have a stable equilibrium point with three species present if a + P < 2 
(ct,P > 0). For the region a > 1 and P > 1 the equilibrium points for all three single species 
are each stable. For the remaining points (a, p) of the parametric space asymptotically 
stable equilibrium points do not exist. 

In (|2.2|l the singularity analysis demonstrates that there is no possibility of an analytic 
solution unless a = P for otherwise r 2 ± are complex. When a = P, r 2 ± coalesce into 



a 



1 



?'2 



l + 2a 

with r 2 a positive integer, n, if 
n + 1 



a 



In - 1 



which is necessarily negative and so beyond the acceptable parameter range of the model. 
Only in the case that a (= P) = 1, for which r 2 = (2) can we expect an analytic solution. 
Then ED is 



(2.3) 



" 1 


1 


1 " 




" M ' 




' " 




" M ' 




l " 




" 


1 


1 


1 




N 







=> 


N 


= h 





+ k 2 


1 


1 


1 


1 




S 









S 




-l 




-1 



where k\ and k 2 are arbitrary parameters, and we have the two constants of integration 
entering at the leading-order behaviour. The linear terms omitted in the analysis of the 
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dominant terms do not cause an inconsistency since they do not affect the leading order 
term of r -2 . 

That the solution is analytic may be demonstrated by the explicit integration of system 
(|1.2|) with a = (3 = 1. When the composed system is integrated and this is substituted 
into the equations for x and y, say, ie z is eliminated using the invariant, the equations 
for x and y decouple and one obtains the solution by simple quadratures to be 



fcie* 
C + e* 



V 



k 2 e t 



1 + (1 - h - k 2 ) e l 
C + e* 



as is suggested by the solutions given in (j2.3j) . 

In terms of the singularity analysis system (|1.2|) is integrable in terms of analytic 
functions at the specific point (1,1) in the (a, (3) plane. This is the point of contact 
between the two regions of stable equilibria reported by May and Leonard ^H] [Fig 1] • We 
observe that the only pattern of leading-order behaviour compatible with the standard 
method of singularity analysis as found in, say, |22| \TA\ is that the exponents of the 
leading-order terms be at —1. However, we may depart 3 from that standard analysis and 
investigate the consequences. If we suppose that p = q = — 1 and r = 0, the dominant 
terms of system (|1.2|) become 



-At 
-Bt 



At" 1 
Br' 1 

c 



(-At- 1 
{(3AT- 1 
aAr^ 1 - 



-aBr- 1 ) 
- Br- 1 ) , 
PBt' 1 ) . 



(2.4) 



The first two of (f2~4")l give 



1 a 

P 1 





A 




1 




A 




B 




1 




B 



1 



1 - a/3 



a 
P 



The third of ()2.4|) gives either C = or aA + (3B = 0. This second condition, coupled 
with the first and second of ()2.4|) demands either that a/3 = l=>a = l,/3 = l, or places 
a and /3 on the circle 



a 



P 



in the (a, P) plane. The former condition, C = 0, coincides with one of the equilibrium 
points with just two species, x and y, present. The other two possibilities, videlicet p = 0, 
q = r = — 1 and p = — 1, q = and r = — 1, correspond to the equilibrium points with two 
solutions given by (0, 1 — a, 1 — P) / (1 — aP) and (1 — P, 1 — a, 0) / (1 — a/3) respectively. 

In a similar situation, if we take the leading-order exponents to be p = — 1, q = r = 0, 
we obtain the leading-order behaviour, 



-At- 2 




At- 1 ( 
B(I3A 



-At- 1 ) , 



3 We do not claim any originality in making a departure. Daniel et al [H] did the same in their study 
of the Heisenberg spin chain with anisotropy and transverse field. For a deep study from the viewpoint of 
cosmological interests see the more recent work by Cotsakis 
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= C (-aAr^ 1 ) , 

for which the solution is obviously (1,0,0). The other two possibilities are (0,1,0) and 
(0,0, 1) , it we recover the equilibrium points with just a single species present. 
Moreover, if one makes a formal expansion 



CiT 




and substitutes this into p.2|> . one obtains (courtesy of Mathematica) that 6j = 0, q = 0, 
Vi, all odd coefficients, a2i+i, are zero and that 

1111 1 1 



GO = 1) a l = ~j fl 2 = TTT' °4 = — T^TT) a 6 = &8 = ~ ^ „„„„„„ , «10 



2' 12' 720 30240' 1209600' 47900160 

which is in accordance with the solution of one species present, videlicet 

' y = 0, z = 0. 



1 — exp[i — in] 

Here we must emphasise again that we are not applying singularity analysis in the sense 
of the Painleve Taste. Nevertheless we see an interesting connection between the results 
of dynamical systems analysis and the simple series substitution. The solution obtained is 
consistent with the Painleve analysis in that it leads to a subsidiary solution j2J although 
the route to its obtention is formally different. Nevertheless it cannot be regarded as a 
subset of the Painleve analysis since a fundamental feature of the analysis is that the 
coefficients of the leading order terms be nonzero. 

We emphasise that this last part of the analysis is not in accordance with the norms of 
singularity analysis as presented in the standard references. Once we admit the possibility 
of a zero as the exponent of the leading-order behaviour of one or more species, we depart 
from the criteria for the application of the Painleve test. Nevertheless we see that results 
can be obtained which are very suggestive and which connect in a natural way with the 
analysis of system (jl.2j) via dynamical systems. For an investigation of the presence of 
two competing species we refer the reader to jllj . 

We conclude our singularity analysis of system Q1.2j) with the final observation that 
generically (|1.2JI is not integrable in terms of analytic functions. Nevertheless the analysis 
has revealed aspects of the properties of the system which perhaps would not be anticipated 
a priori. 



3 Symmetry analysis 

The system (|1.2|) is autonomous and so possesses the Lie point symmetry dt- For integra- 
bility in the sense of Lie we require the knowledge of a three-dimensional solvable algebra. 
The knowledge that (jl.2ft is a system of first-order differential equations and so possesses 
an infinite number of Lie point symmetries does not help us to find the additional two 
symmetries. For the purposes of symmetry analysis we make a change of variables to 
convert system IJ1.2JI to a quadratic system 4 . 

4 For a system of the general form of 11.11 this is not possible, but the simplifying assumptions of May 
and Leonard that = r — * 1, i = 1, n, under a rescaling of time does enable the transformation of system 
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We write 



X = xe-\ Y = ye- 1 , Z = ze~\ T = e*. 

Then system (jl.2j) becomes 

X' = -X (X + aY + f3Z) , 
Y' = -Y (f3X + Y + aZ) , 
Z' = -Z (aX + /3Y + Z), 



(3.1) 



where we use the prime to denote differentiation with respect to the 'new time', T. If one 
assumes a solution of the form X oc T p , Y oc T q and Z oc T r , one finds that 



X(T) 
Y{T) 
Z{T) 



1 + Q + /3 



(3.2) 



is a solution. It is not the general solution because it does not depend upon three arbitrary 
constants. It is simply a peculiar solution. In fact, going back to the original variables 

X = xe~\ Y = ye~\ Z = ze'\ T = e*, 

solution (|3.2|) corresponds to the equilibrium solution 

(1,1, 1) T 

of the original system (|1.2|) . Note also that to every equilibrium solution of the original 
system (|1.2|) there corresponds a peculiar solution of (|3.2jl . 
By inspection (|3.1jl possesses the two Lie point symmetries 



Ti = dr, r 2 



-Td T + Xd x + Ydy + Zd z 



(3.3) 



with the Lie bracket [ri,r2] iB = —V\. 

In analogy with (|1.4|) we add the constituent equations of ()3.1|) to obtain 



(X + Y + Z)' = - {X 2 + (a + (3)XY + Y 2 + (a + [3)YZ + Z 2 + (a + [3)ZX] . (3.4) 

In the particular case that a + (3 = 2 a possible source of additional symmetry is from the 
decomposition of symmetries of the composed system of ()3.1(l . videlicet 



N' + N 2 = 0, 



(3.5) 



where N = X + Y + Z . We are unaware of this approach being used in the literature 
before now. 



11.211 to the simpler form. The change of variables used here is not beneficial for the singularity analysis 
the outcome of which is very much dependent upon the representation of the coordinates used, but, as it 
is a point transformation, has no effect upon the algebraic structure of the system. 
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As a first-order differential equation (|3.5j) has the same problem of determination of 
symmetries as the composed system (j3.1j) . but it can be written as the second-order 
differential equation 

w" = 0, (3.6) 
by means of the Riccati transformation 

N = ?L w = exp (^J NdT^j , w' = iVexp (J N&T^j . 

A symmetry E = Qdr+^d w of ()3.6j) can be written as a symmetry of (|3.5|) . A = Tdr+Tjdjsf, 
as follows. The first extension of E 

S [l] = QQ T + E Q w + ( H / _ ^/Q/) Q w , QQ T + f£ (^j _ N Qf\ 



Si 


= <9 W 


s 2 


= rd w 


E = 


= wd w 


s 4 


= dr 


E 5 


= 2Td T + wd w 


E 6 


= T 2 d T + Twd w 


Ey 


= wOt 


s 8 


= Twd T + w 2 d w 



so that 

The Lie point symmetries of (|3.6j) transform as follows 

+ Ai = N exp (- f NdT) d N 

* A 2 = (1 - TN) exp (- / AdT) d N 
lost (The source of the Riccati transformation) 

+ A4 = 8t 

* A 5 = Td T - Nd N 

* A 6 = T 2 d T + (1 - 2TN) d N 

* A 7 = exp (/ AdT) (dr - N 2 d N ) 

* A 8 = exp (/ AdT) (Td T - TN 2 d N ) . 

In A4 and A5 we have the T\ and r 2 of (|3.3j) . If we examine the remaining symmetries, 
A 2 and A@ do not decompose. The remaining symmetries decompose according to 

Ai -> A' 11 =exp(-fNdT){Xdx + YdY + Zdz + 2(X'dx>+Y'd Y > + Z'd z ,)} 

A 7 -> A 7 1] =exp(/iVdr){9 T -iV[Xax + ^ + ^ + 2(A / ax' + ^ / 5y/ + Z'a Z /)]} 

A 8 -> Ag 1 ' = exp (J A"dT) {Tdr — TN [Xdx + i^cV + + 2 (X'dx 1 + + Z'dz 1 )]} , 

where we have written the first extensions of the decomposed symmetries to highlight the 
discomforting point that these three symmetries bring no new information. The effects of 
Ai, A7 and A 8 are the same as that of T 2 Ag) on the autonomous system. We conclude 
that it is possible to decompose symmetries just as it is possible to decompose equations, 
but the results are not necessarily useful. 

We already have an invariant derived from the composition of system (jl.2p with a +/? = 
2 in (|1.6|) . videlicet 

t 1 T 

1 X+Y+Z ' 
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when written in the new coordinates. It is evident that I\ is not an invariant associated 
with either T± or T2 and so we may use T± and T2 to seek a new invariant, in fact an 
integral if we require the function to vanish under the action of both symmetries 5 . 

In essence we use the method of reduction of order |171ll8j with the two symmetries T\ 
and IV The former is a consequence of the autonomy of system (|3.1j) and we eliminate T 
as the independent variable in favour of Z by writing 6 

dX _X(X + aY + /3Z) dY _Y (PX + Y + aZ) 

~dZ ~ Z (aX + (3Y + Z)' ~dZ~ Z (aX + /3Y + Z)' ( ' ' 

System (|3.7j) is homogeneous with the obvious symmetry T2 = Xdx+Ydy + Zdz following 
from T2- Under the standard change of variables 

X = uZ, Y = vZ, f 2 -> Zd z 

we may eliminate the ignorable coordinate Z (actually in the form exp (-Z)) to obtain 
the single first-order differential equation 

dv _ v [(/? - a) u + (1 - 0) v + (a - 1)] 

du u [(1 - a) u + {a - 0) v + (fi - 1)] ' 1 ' ' 

It is a trivial matter to integrate Q3.8|) in the case a = (3 which, we recall, is the condition 
for the singularity analysis to give real resonances. The integral is 

u (v — 1) ' 

in which we note that the parameter a is absent. 

When a ^ P, the first-order differential equation is not integrable in closed form (as 
far as Mathematica is concerned). When a + /3 = 2, the integration of (|3.5|) gives 

h = ~ T (3.9) 

which corresponds to I\. 

May and Leonard note that the product xyz (in our notation) has an interesting asymp- 
totic behaviour. 

From (|3,lj) we find that 

(XYZ)' = -XYZ{l + a + P){X + Y + Z) (3.10) 

so that 

(XYZ)' 



XYZ 



-kN 



5 Although one usually looks for an integral/invariant associated with a a single symmetry - the only 
way possible for a two dimensional system - there are at times great benefit and simplification to imposing 
the requirement that the integral/invariant be associated with two (or more) symmetries [1]. 

6 There is no essential difference made by the particular choice of a new independent variable. 
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and in the case that a + (3 = 2 so that 1)3.5)) and 1)3.9)1 apply this is 
{XYZ)' _ N' 

XYZ ~ 3 77 

whence 

XYZ 

h (3.11) 



(X + Y + Z) 3 

which also may be written as 

XYZ = I 2 {Ji + T) 3 . (3.12) 

The form (|3.11|) indicates that this integral corresponds to the symmetry dr- The 
group theoretic origin of 1)3. 11 j) is easily seen. The invariants of dr are X, Y and Z. 
The requirement that f(X,YZ) be an integral of ()3.1|) leads to the associated Lagrange's 
system 

x 7 ~ Y 7 ~ ~z 7 ' 

We may use the theory of first-order differential equations [H] [p 45] to combine the elements 
in a specific fashion to give, when ()3.1j) with a + (3 = 2 is taken into account, 

dX dY dZ d(XYZ) 



X' Y' Z> -3XYZ(X + Y + Z)' 
Taking, say, the first with the fourth and using the equivalence of dX/X' to dT we have 
d(XYZ) 



XYZ 



-3NdT 



and ()3.11|) follows when the composed equation, 1)3.5)1 . is used. 

We may use ()3.9)l and ()3.12)) to eliminate Y and Z (say) from system ()3.1)) with a + (3 = 
2. To maintain a certain compactness of notation we write X + Y + Z as N(T) and XYZ 
as m(T). We obtain 

Z = ^- and Y 2 + (X-N)Y + - = (3.13) 
X i X 

so that 



Y = \^N -X±^{N - Xy (3.14) 
The first-order differential equation satisfied by X is found to be 

X' = -aX T |(a - /S)xJ (a - X) 2 - ^. (3.15) 

The single first-order differential equation, 1)3.15)1 . for X(T) encapsulates the informa- 
tion already gleaned by our various analyses. The contents of the root ensure our inability 
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to find a solution to 1)3.15)) unless a = (3. When a = [3 = I, 1)3.15)1 is always integrable 
since then (|3.15|) is simply a linear first-order equation. 

It would be evident to a reader with only a modest acquaintance with symmetry analysis 
that the procedure of this Section is not unique. We could look to replace the system 1)3.1(1 
by either a single third-order differential equation or a second-order differential equation 
plus first-order differential equation which is a standard part of the method of reduction 
of order [12 HE]- 

For the system ()3.1() the former option is not feasible as leads to a very awkward 
algebraic equation. The latter option is feasible and is found in the analysis of the Euler- 
Poinsot system )19) . Indeed it is a logical consequence of the case for which two equa- 
tions are replaced by a single second-order ordinary differential equation )12) 114) 120) . In 
the papers cited the symmetries of the second-ordered differential equations sought were 
Lie point symmetries. In their study of the two-dimensional predator-prey system with 
malthusian growth Baumann and Freyberger [2] replace their two first-order differential 
equations with a polynomial second-order differential equation and then seek generalised 
symmetries with specific structure of the second-order equation. These give generalised - 
equally point since the two coalesce for systems of first-order equations - symmetries of 
the original Lotka-Volterra system. From the symmetries integrals and invariants follow 
easily. 

We have not follow the procedure of Baumann and Freyberger in this paper since our 
investigation is, as the title of the paper proclaims, of analytic solutions and the range 
of parameters is already set by the singularity analysis. Were our intentions otherwise, 
a symmetry analysis along the lines of those in the papers cited above would be appro- 
priate. The problem with (systems of) first-order equations is that the number of Lie 
point symmetries (equivalent, as noted above, to generalised symmetries) is infinite and 
so there is no finite algorithm for their determination. The increase of order, which is an 
integral component of the method of reduction of order, makes it possible to implement 
a finite algorithm. However, if like Baumann and Freyberger one introduces generalised 
symmetries at the higher order, the finite algorithm is lost. One may as well substitute 
Ansatze of choice into the original system of first-order equations. 

Our concern with the analytic behaviour of the system led us to use a rather restricted 
symmetry approach to the determination of appropriate symmetries. If one removes the 
requirement of analycity, a wider investigation, even as general as that of Baumann and 
Freyberger, would be appropriate. 

4 Conclusion 

Our investigation of the model for competition among three species presented by May and 
Leonard was motivated by the singularity and symmetry analyses which are appropriate 
to integrable systems. We found that the critical values of the parameters revealed in the 
analysis of the system using the methods of dynamical systems were echoed in subsequent 
analyses from the viewpoint of the singularity and symmetry approaches. That (1), asq 
()3.1)) . is a decomposed system when a + (3 = 2 made our analysis easier since the composed 
equation, (|3.5j) . is trivially integrable. This provided one invariant for system ()3.1I) . The 
second invariant was suggested by the analysis of May and Leonard who showed that 
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XYZ was expressible in terms of an invariant for all values of a and [3 provided that the 
populations were small. When a = (3 = 1, the asymptotic invariant becomes a global 
invariant. The existence of these two invariants, corresponding to the two symmetries 
of invariance under time translation and similarity transformation (more evident in ()3.1|) 
than (1)), enabled the system of three autonomous first-order differential equations to be 
reduced to a single nonautonomous first-order differential equation. The consequence of a 
lack of further symmetry for general values of a and (3 even with the constraint a + (3 = 2 is 
quite evident in the form of the equation with a + (3 = 2. The system (|3.1j) cannot exhibit 
chaos since the autonomous invariant 1)3.11(1 implies that the system can be reduced to 
an autonomous system of order two and is thereby formally integrable. However, this 
integrability is not in terms of an analytic function let alone in closed form 7 . This is quite 
evident from the form of 1)3.15(1 . The improvement in the integrability of 1)3.15(1 as the 
parameters become more closely aligned to the values which give favorable results for the 
singularity analysis is clearly apparent. The simple removal of the imaginary part of the 
resonances when a = f3 enables ((3.15J) to be integrated trivially in terms of an analytic 
function. In this respect the system 1)3. ljl is an excellent paradigm for the implications of 
the requirements of the singularity analysis. 

May and Leonard rightly indicate the marked change in the behaviour of the system 
as the dimensionality is increased from two to three. The numerical results which they 
present for a + (3 > 2 are very suggestive of the behaviour in the solution which one would 
expect when the resonances are complex. Indeed one would expect similar behaviour for 
a + (3 < 2, a ^ (3, but the change in sign of the real part of the exponent means that 
there is damping of the oscillations rather than growth. Indeed, if one considers the basis 
of the model, the values of the parameters a and (3 should be such that a + (3 is likely to 
be an order of magnitude less than two. 

The community matrix of system (l)/ l)3.1j) is circulant and a general expression for 
the eigenvalues is available. When a = f3, the community matrix becomes symmetric. 
Although May and Leonard make the point that one could scarcely be interested at looking 
at the analysis of a system of dimension greater than five, there may be some merit in 
the consideration of systems of dimension greater than three. If all of the nondiagonal 
elements are equal, which is a fairly drastic extension of the a = (3 case discussed above, 
all submatrices containing the principal diagonal are circulant and the possibility of the 
existence of composed systems leading to invariants is real and integrability for competition 
among n species is conceivable. The constraint on the interaction coefficients to make the 
community matrix circulant and so system (1) amenable to some analysis may be regarded 
as severe. However, as we observed in the Introduction, among species of similar habit 
as well as habitat the interaction coefficients are likely to be less dominant than the self- 
specific effects. This is not a small step from assuming equality. 



7 We do not into a discussion of the meaning of integrability in terms of functions which are not analytic. 
Although the formal requirement that a solution be analytic is accepted, there are sufficient acceptable 
exceptions for a certain laxness in practice. 
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